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I. " INTRODUCTION 


The purpose of this investigation was to develop a heat 
transfer model that will illow prediction of the temperature 
distribution in a container stored rocket motor placed in a 
hostile thermal environment such as the desert. It is pro-~ 
pesed that such a model would be a useful tool for thermally 
optimizing future contaiier designs. As extreme variations 
in the rocket motor temperature may lead to large thermal 
stresses in the propellant which could result in fracture, 
or otherwise degrade the performance of the motor, a major 
objective of this study was to design a model that could 
reliably predict the thermal gradient in the motor. The 
predictions would be based on the surface temperature dis- 
tribution, the thermal properties and the geometrical de- 
tails of the system. The model aay also bevusnd eG predict 
a critical temperature range over which the propellant must 
be chemically stable when in a storage situation. The 


upper limit of this temperature range is referred to as the 
design temperature of the system. As the design temperature 
for most weapon development projects is derived from dump 
storage eonditicns; a dump storage situation was used to ob-~ 
tain the experimental data for this project. 

Several approaches were taken to predict the rocket 


motor temperature distribution from a knowledge of only the 


surface temperature distribution of the storage container 


and the thermal properties and geometrical details of the 
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experimental model. The experimental model used in this 
test was a once-fired Navy antisubmarine rocket (ASROC) 
- motor, filled with dry desert blow sand to simulate the pro- 


pellant, and placed in its storage container. This container 


opt 


system was placed in a dump storage site at the Naval Meanone 
Center, China Lake, California to simulate a desert environ- 
ment. 

The method of complex temperatures [Ref. 1] was used to 
develop an analytical prediction of the transient tempera- 


ture field that exists in a container stored rocket motor. 
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The analytical model assumes that heat is transferred only 


ea aad 


{ in the radial direction and that the container surface tem- 


perature variation is sinusoidal with time. Comparison 


between theory and experiment is within experimental uncer- 
tainty when temperature is interpreted as "bulk" temperature. 
The analytical model is especially useful for studying 
géometrical and thermal physical property effects on rocket 
motor temperature. Such pasvameter studies have been carried 
out and the results are presented in a form that will be 
isetal fron a container design point of view. 

TRUMP [Refs. 2 and 3], a computer program for transient 
‘and steady-state temperature distributions in multidimensional 
| systems, was used to obtain detailed information about the 
| : thermal state of the rocket motor. TRUMP allows actual con- 
- tainer surface temperature distributions to be used as well 
as sinusoidal variations. In addition, both one dimensional 


(radial) and two dimensional (radial and circumferential) 
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heat transfer were modeied with TRUMP, using both the sinu- 
soidal and actual temperature distributions. The actual . 
temperature distributions were obtained from the experimental 
data of the motor container systen, 

Compartsons between the experimental values and those 
predicted by the models were in good agreement ,with those 
predicted by TRUMP using the actual temperature distribution 
as the boundary condition being the closest. However, the 
sinusoidal variations used in both the analytical model and 
the TRUMP model are also suitable for design purposes. 

Another aspect of this project was to obtain the storage 
container surface temperature distribution using cholesteric 
liquid crystals, a material that undergoes brilliant changes 
in color over known, well defined temperature ranges. Color 
slides and movies were taken of the liquid crystals demon- 
Strating the feasibility of using them for on site tempera- 


ture measurements. 
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LI... BACKGROUND 


In 1959 the Naval Weapons Center, China Lake recognized 
t the need for a concerted attack on the problem of chexwal 
criteria assignment for new weapon systems. In 1963 a task 
force was established to study the complete environmental 
criteria determination problem. The key to this problem 
seemed to be the thermal area in the storage and transporta- 
tion events of any item. It was realized that transportation 


was a short term situation cempared to the storage situation. 


in storage. There are three types of storage; covered, igloo 
and dump. The dump storage situation leads to the more 

£ extreme thermal exposure situations which then leads to the 
design temperature. 

As data was not available for the dump storage situation, 
instrumented storage dumps were created at representative 
piaces ona woridwide basis so that statistical data could 
be derived on a variety of ordnance. The first site was at 
China cakes California, in the middle of the Mojave Desert. 
This site now has the capability to return about 250 channels 
of information on a continuous time~temperature basis (Fig- 


ure 1). Other arctic and tropical sites were set up to study 


Therefore, the major portion of the life of an item must be | 
| 
! 
t 
i 
: 
i 


extreme conditions. 
The dump storage situation was reproduced to study the 
extreme situatiion. The ordnance was exposed singly, directly 


situated on the ground, with the long axis aligned in the . 
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north-south direction to allow maximum normal exposure of the 
container surface to the sun's rays. In actual practice, 
ordnance is usually stacked and oriented in other than a 
north-south direction, thereby svoldine the extreme situa- 
tion. Ordnance sitting on the ground receives reflected 
radiation from the ground, cannot quickly give off heat by 
conduction to the soil, and its not as apt to be cooled by 

the prevailing breeze; therefore, extreme temperatures 
result. 

The most important source of heat to the ordnance is the 
direct radiation from the sun, with reflected radiation of 
secondary importauce. For extreme conditions to occur the 
wind must be caim (less than 5 knots), the sky clear, and 
the outside air temperature high. After sunrise, the ord- 
nance skin temperature rises much more rapidly than the 
ambient air temperature; therefore, the surrounding air 
cools the ordnance, rather than heats it. 

The rocket motors used for the tests were military sur~- 
plus. Even though the material had served its intended in- 
Fleet purpose, it was still deneScaneative of new -hardware, 
when teed in a thermodynamic context. When inert rocket 
motors were available, they were used intact; however, in 
most cases, once-fired hardware was used. Thoroughly dried 
desert blown sand, being similar in thermal properties to 
most propellants, was used to backfill empty rockei motors. 
It was assumed that the thermal response of the sand filled 
motors was essentially the seus as actual propellant filled 


motors. 
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IIL. EXPERIMENTAL PROCEDURE 


Although Naval Weapons Center, China Lake had accumulated 
vast amounts of data a the past, it was decided to in ‘u- 
ment a rocket motor storage container system especially fox 
this project. This would allow base data to be taken exactly 
where it was required. It also allowed variations in the 
system without interfering with one of China Lake's ongoing 
projects, An ASROC system was chosen for this study. The 
outer storage container was 75 inches long with an inner 
diameter of 18 inches and a wall thickness of 1/16 inch. The 
rocket motor was 57 inches long with an outside diameter of 
12 inches and a wall thickness of 1/4 inch. Both the con- 


tainer and motor were made of steel. 


The rocket motor storage container system was instrumented 


with 20 gage copper-constantan insulated thermocouple wire 
which has an ISA Calibration of + 1-1/2°F over the range ~75 
to +200°F. Twenty-one thermocouples were originally placed 
on the syebouswith-post tions indicated in Figures 2 and 3. 
The ambient air temperature was measured with thermocouple 
number 19 which was located in a Stevenson shelter about 60 
feet away from the system (Figure 4):. The thermocouples were 


mounted intrinsically on the motor and storage container. 


. Two small holes were drilled approximately 1/8 inch apart in 


the metal and the individual wires were inserted in the holes. 
The metal was then hammered around the wires until a snug 


fit was obtained. Bead thermocouples were mounted at the 
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Figure 2. 
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Thermocouple Locations on Experimental Systen. 


thermocouples were located under the section painted 


Their locations corresponding to 


shown above are: #14= #1, pee #2, #16= #8, #17= #9, 
and #18= #3 (See Figure 3). 
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Figure 4. Stevenson Shelter. 


ee. tas 


peer sire 


center of the motor and in the air gap. The thermocouples 
located at the center of the motor were supported by small 
pieces. of wood several inches from the head. The use of 
these supports was necessary to keep. the thermocouples in 
position when the motor was being filled with sand. After 
all the thernocouples on the rocket motor were in place, 
the rocket motor was filled with dry desert blown sand. 

The wires from the two thermocouples located in the center 
of the motor were led out a hole in the end cap. To avoid 
settling of the sand after the motor was in place on the 
site, with a resulting air gap being formed between the 
sand and the motor skin, the sand was compacted by striking 
the sides of the motor with small sledge hammers and then 
adding additional sand through the hole in the end cap. 
This was continued until the sand was tightly packed. The 
hole in the end cap was then sealed. The rocket motor was 
carefully placed in its storage container (Figure 5) which 
had previcusly been instrumented with thermocouples. The 
thermocouples in the air gap were mounted by affixing the 
lead wire to the rocket motor at the desired position and 
then putting a 90 degree bend in the wire so that it placed 
the bead of the thermocouple approximately 1.5 inches into 
the air gap. Neither the thermocouples in the center of 
the motor nor those in the air gap could be considered accu- 
rately positioned; however, every effort was made to mini- 
mize Seeteaenine errors. All thermocouple wires were 


located inside the storage container and were ted through a 
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; hole in one end. This hole was then sealed. The two halves 
of the storage container were then bolted shut. 

The outer surface of the rocket motor and the inner and 
outer surfaces of the storage container were all painted 
various shades of haze gray. Weathering had caused the 
painted surfaces to appear fairly rough. This is typical 
of the conditions of a storage dump. From the condition of 
the surfaces, it was estimated that the emissivity was 
approximately 0.9. 

Prior to loading the rocket motor into the storage con- 


tainer, it was decided to apply liquid crystals (See Appen~ 


dix A) to part of the storage container surface in order to 
obtain a thermal mapping of the surface temperature at any 
instant of time. Liquid crystals are temperature sensitive 


materials. that produce immediate thermal images in a pattern 


of colors which respond rapidly to minute changes in sub- 
strate surface temperatures. A second reason for applying 

Hy the crystals to the container surface was to determine the 
feasibility of using the crystals under adverse environmen- 
tal conditions (desert dtaosphere): Prior to applying the 
scgevaies a 15 inch strip of the storage Scntainer; 20 inches 
from one end, was sprayed with two coats of Testors Spray 
Pla E..amel No. 1249, Flat Black as a background for the 


crystals. A one inch strip of 11 different ranges of crys- 


tal, with approximately 1/2 inch of black paint between 


ere: 


them, was applied over the black paint. Two coats of each 
crystal were applied, using a small paint brush. The first 


coat was allowed to dry completely before the second coat 
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was applied. After the crystals were dry, two coats of Rez 
polyurethane (gloss clear plastic coating, interior-exterior 
77-5) coating were applied by. brush completely sveehe the 
crystals and black urea area. The polyurethane coating 
was applied to protect the crystals from wind blown sand and 
from the ultraviolet rays of the sun. Ten of the eleven 
crystals had been previously calibrated [Ref. 4]. Using the 
constant temperature bath procedure recommended in Ref. 4, 
R-27 was calibrated and the compli :e calibration results are 


shown in Table I. 


The rocket motor storage container system was then moved: 


to the China Lake dump storage site. The system was aligned 


in a north-south direction, well away from the influence of 
other ordnance (Figure 6). The thermocouple leads were 
connected through a junction box and underground cable to a 
Honeywell Electronik 25 Recorder which had been calibrated 
to read the thermocouple output aivectiy in degrees Fahren- 
heit to an accuracy of ae 1°F. The recorder was located in 
an air-conditioned building about 60 feet from the system. 

Initial data indicated that the number 7 thermocouple 
was not responding properly and therefore this data was 
neglected. Initial color photographs were taken of the 
liquid crystals and it was immediately apparent that good 
thermal mappings could be obtained if the crystals were 
stable under the adverse desert environment. The brilliance 
of the colors exhibited by the crystals under the bright 


desert sun was much better than had been expected. The 
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TABLE L 
Calibration of Liquid Crystals 


Calibration Bath 


NCR Manufacturer's 2 Coats 
“Desig. ‘Color Change —.__—Responses ss Liquid ‘Crystals 
; *C °¢ 
Red 27.0 25.6+.5 
Green 28.6 28.0+.5 
Blue 30.0 28.7+.5 
R-33 
Red 33.0 32.7+.5 
Green 34.6 33.3+.5 
Blue 36.0 34.2+.5 
R-37 
Red 37.0 36.24+.5 
Green 38.6 37.1+.5 
Blue 40.0 38,0+.5 
R-41 . 
Red 41.0 40.345 
Green 42.6 41.0+.5 
Blue 44.0 42.0+.5 
R-45 
Red. 45.0 42.8+.5 
Green 46.6 43.6+.5 
R~-49 | 
Red 49.0 46.7+.5. | 
Green 50.6 47.1455 ; 4 
Blue 52.0 48.44.5 i 
R-53 ] 
Red 53.0 50.54.5 4 
Green 54.6 52,.14+.5 j 
"Blue. 56.0 53.34.5 4 
R-56 j 
Red: 56.0 53.84+.5 i 
Green. 57.6 56.0+.5 
Blue 59.0 56.5+.5 
R-59 | 
Red 59.0 56.9+.5 j 
Green 60.6 57.5+.5 
Blue 62.0 58.94+.5 
S-62 e 
Red: 62.0 60.1+.5 
Green 62.6 : 60.4+..5 
Blue 63.0 60.9+.5 
S-64 4 
Red 64.0 60.9+.5 ; 
Green 64.6 61.4+.5 4 
Blue - 65.0 62.7+.5 : 
a 
4 
25 : 
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Experimental System at Dump Storage Site. 


Figure 6. 
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system was allowed two weeks to reach a periodic steady 
state before additional photographic data was obtained. 

Extensive photographic data was collected on 27 and 28 
July 1972 after two weeks of exposure to the desert environ- 
ment. Both super 8 mm and 16 mm color movies and 35 mm color 
slides were taken of the liquid crystals. No colored filters 
were used on any of the cameras, although standard haze fil- 
ters were used to take the super 8 mm movies and most of 
the 35 mm slides. 

At this time, a second storage container, this one with- 
out a rocket motor inside, was instrumented with intrinsic 
thermocouples in the same manner as-the previous container. 
As only three data channels remained open on the recorder, 
only three thermocouples were applied to this new container. 
The three thermocouples were applied at the 0300, 0900, and 


1200 positions at the midpoint of the container. This con- 


tainer was set end to end with the system that was already 

in place at the site. The purpose of this study was to 
determine if the inclusion of the rocket motor in the con- 
eines had a significant effect on the surface temperature 

of the container. Thermocouple #7 was connected at the 0900 
position, #23 at the 1200 position, and #24 at the 0300 posi- 
tion. It was immediately apparent that thermocouple #7 was 
continuing to give unreliable readings and therefore the 


data taken on channel #7 was neglected. 
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IV. THEORETICAL ANALYSIS 


A. ONE-DIMENSIONAL ANALYTICAL MODEL 
The first step was to try to devise an analytical model 
that would simulate the actual rocket motor storage contain- 
er experimental system. The first simplifying assumption | 
was that the storage container temperature could be modeled | 
by a sine wave which had a period of 24 hours. A compari- 
son of the sinusoidal variation to the average (bulk) stor- 
age container temperature [obtained by averaging the four | 
thermocouple readings on the surface of the container (1, 
8, 10, and 12) as shown in Appendix D] is given in Figure 7. 
The method of complex temperature as presented by Arpaci 


{Ref. 1] was used to find the steady periodic solution of a 
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body experiencing a periodic sinusoidal disturbance. A com- 


plete analytical derivation is given in Appendix B. The 
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gcneral heat conduction equation in cylindrical coordinates 


was the basis for this derivation. It was assumed that 


there was one dimensional radial heat flow with no conduction 


in the axial or circumferential directions, that no heat 


sources existed in the model, that the rocket motor storage 
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container system was infinitely long, and that the sinusoid- 


ne 


ally varying surface temperature was spacially uniform over 


ere nce 


the entire container surface. The storage container tempera- 
ture is assumed to vary as 


T, = (y- T,) sin wt + 3, 
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where T,, = maximum bulk temperature of the storage 

container 

T, = average bulk temperature of the storage 

container 

® = frequency of tne sinusoidal variation 
(214/24 hours) 

t = time (hours) 


It was assumed that all the thermal properties remained 
constant over the temperature range of the problem. The 
effective heat transfer coefficient, h, across the air gap 


between the storage container and the rocket motor combines 


duction into one coefficient. The radiation coefficient was 
linearized by assuming constant representative temperatures 


in the equation 


_ 2 2 
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RAD 2 


where o is the Stefan-Boltzmann constant and is the 


1-2 
radiation exchange factor. The convection coaEFicient is 
the effective conductivity of air, obtained from the Beck-~ 
mann correlations [Ref. 5], divided by the width of the air 
gap. In the analytical model, the effective conductivity 
was assumed to equal the conductivity, thereby treating it 
as pure conduction and giving the equation 

heh, 
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the heat transfer effects of radiation, convection, and con- | 
| 
{ 
| 
| 
N- | 


An initial condition was not specified in this deriva- 
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k 
tion as the only concern was the steady-state, periodic be- | 


havior. The steady~state solution is (Appendix B) 


T(x,t)-T,  YBER,* (af) +BE4,,? (aE) 
6(r,t) = =P sd poi 
M” A Vx? + X,? 


R i 


sin(Wt+é*) 


O(r,t) = 6 sin (wt + 6*) 
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where T(r,t) = the temperature of a point r in the rocket 
motor at time t 


wr ? 


= conduction parameter 


Ee = = dimensionless distance from the center 
o of the rocket motor 


rot inner radius of the rocket motor 
r = distance from the center of the rocket motor 
a= se = thermal diffusivity 
p = density 
k = thermal conductivity 
C = specific hest 
BER = real Bessel Function 
BEL = imaginary Bessel Function 
X, = BER (a) + —*— BER, (a) + —*— BEi., (a) 
R re) /28 1 /2B 1 
X, = BEi_ (a) + —\— BEi,(a) - —* — BER, (a) 
i re) /2B 1 /2B 1 


ir, 
6 = ra Biot modulus 


1 BEA, (a&)X, - BER, (ab)X, 


6* = tan ee 
BER (ab) xX, + BEL (ab) X, 


Two computer studies were done based on the steady state 
solution. The first study was a completely dimensionless 
situation which served as a parameter study of the effects 
of varying a and 8 on the temperature and the time lag of 


the temperature at various positions in the model. 


- 


= conduction parameter 


and . 
hr, 
8 =—— = Biot modulus 
Parameter a was varied from 1.0 to 5.0 and 8 was varied 


from 0.1 to 100. These were the only values studied, as 
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only values within this range are of interest in this type 
problem, The computer program and its output are given at 
the end of Appendix 5. The output lists the following 
values: 

1) a, the conduction parameter 

2) 86, the Biot modulus 


3) &, the non-dimensional distance from the center of 
the motor 


4) 6%,the time delay between the maximum storage con- 
tainer temperature and the maximum temperature 
reached at the point of interest in the motor 

5) 6_,the relative amplitude of the maximum tempera- 
ture at the point of interest compared to the 
maximum temperature of the storage container 

The time delay is given in radians, where 2m radians equals 
one complete cycle. A graph of the time delay versus 8 for 
a constant value of “a" is given in Figure 8 at three dif- 
ferent positions within the motor. A graph of the relative 
amplitudes of the temperatures versus 8 for a constant value 
of "a" igs given in Figure 9. It was noted that for a con- 
Stant value of "a", the time delay decreased as 8 became 
larger. As the point of interest approaches the center of 
the rocket motor, the time delay increases. The relative 
amplitude of the temperatures also becomes larger as 8 is 
increased when the value of "a" is held constant. If 8 is 
held constant and "a" is varied, the time delay increases 
and the relative amplitude decreases as "a" increases. 


The second study was obtaining the analytical solution 


to the particular rocket motor storage container system. 
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studied at China Lake. The thermodynamic properties of dry 
sand were obtained from Ref. 6 as 
94.8 lbm/£t> 


0.188 BTU/hr ft °F 
0.195 BTU/1bm °F 


p 
k 
¢ c 


Substituting these values and using 1440 minutes (24 hours) 
as a complete cycle, the parameters a and £ were calcu- 
lated for this model as 


oO 
ae 7 2.43 


where bai 5.75 inches, the inner radius of the rocket motor, 


Er eee 


where h = h,_.. +h 


con * PRap 
k 

and hoon = woe e626 x 10°" nt 
3 hr-f£t°°F 


where Ar = 2.94 inches, the distance across the air gap 


~2 _ BTU: 
and k,jp = 1.62 x 10° (Fase 
2m 2 BIU 
hoan = y_50(Ty+ Ty) (T,74T,“)= 1.09 
RAD pe eles Naat de aes a3 hr £t2°F 


where o is the Stefan-Boltzmann constant, F is the 


1-2 
radiation exchange factor which for this gecmetry is 


1 
ca eae Ried. eee = 0.84 
ce, ce, 
L742 


when €, = & = 09, x," 6.0, ro * 8.94 


: therefore fh = 1.15 BTU/hr £t2°F 
The average surface temperature of the storage container 


was found to be 104°F for a particular day at China Lake, 


with a maximum temperature of 138°F, These values were ob- 


tained by averaging the readings of thermocouples 1, 8, 10, 


ture. 

The temperatures of seven positions within the rocket 
motor were calculated and the results are printed at 30 
minute intervals for one complete cycle in Appendix B. A 
graph of temperature versus time was plotted by the computer 
showing the relationship between the surface temperature of 
the storage container (TINF), the temperature on the outer | 
skin of the rocket motor (TEDG)-., and the temperature at the i 


center of the motor (TCEN). This graph is Figure 10. 


and 12 as shown in Appendix D which give the bulk tempera- | 
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B. TRUMP MODEL | 
The rocket motor storage container system at China Lake 
was modeled on TRUMP, a numerical conduction code, (See 
Appendix © for a description of the TRUMP program) to pre- 
dict the temperature at any point in the system from a know- 
Tedge of the storage container surface temperature variation, j 
the thermal properties and the geometrical details of the 
{ 
system. Two models were used to simulate the rocket motor 
storage container system and several variations of each 
model were investigated. 


The first model assumed one dimensional heat transfer 


(radial). The system was modeled as two infinitely long 
concentric steel cylinders, the inner of which was filled 
with dry sand. A 2.94 inch air gap separated the cvlinders. 


The model was subdivided into concentric volumetric elements 


| 
| 
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with representative nodal points as given in Figure 27, 
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Appendix C. It was assumed that the storage container sur- 


face temperature was spacially uniform. From the data given 
‘ in Appendix D and the observation of the liquid crystals' 
thermal mapping, it was obvious that the temperature distri- 


bution on the storage container was not spacially uniform. 


In order to simulate a spacially uniform condition, the 


readings of the thermocouples located at the 1200, 0300, 


and this average value of the surface temperature (referred 


we o 
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to as the bulk temperature) was used as the spacially uniform 


0600, and 0900 positions (#1, 8, 10, and 12) were averaged 
' 
{ 


temperature distribution. Two methods were used to describe 
the container temperature. The first method: used the maxi- 


mum bulk temperature (138°F) and the average bulk tempera- 
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ture (104°F) of the storage container to generate a sine 


wave with a period of 24 hours (1440 minutes). The second 


vals and fed this data into the TRUMP program in a tabular 
(temperature versus time) form. The version of TRUMP used: 


| method took the bulk temperature readings at two hour inter- 
in this problem was limited to a table length of 12 tabular 
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values. TRUMP interpolated between the tabular points. 
Figure 11 compares the actual bulk data with the sinusoidal 
approximation and the interpolated tabular values. 


Several assumptions were made to simplify the solution 
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F of this problem. As the thermocouple data from the storage 
container gave an average value of the temperature across 


the 1/16 inch steel wall, node 12 was modeled as a zero 
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volume boundary node with a known temperature impressed upon 
it. It was also assumed that heat transfer across the air 
gap occurred by radiation and conduction alone. Free convec~ 
tion effects were initially neglected. This assumption was 
later modified to investigate the free convection effects. 
All surfaces of the storage container and the outside surface 


of the rocket motor were painted various shades of haze gray 


id arte ta ARTE A 
*- unio 1 CoM “ 
+ 


and it was estimated that the emissivity of these surfaces 
was 0.9. The radiation exchange factor, F ug ttor this 


model was the same as that for the analytical solution CFs 


y 
3 
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= 0.84). It was aie assumed that there was perfect thermal 
contact between the rocket motor and the sand that filled: 
it. This neglects the possibility that the sand might 
slightly settle after being on the site for a long period of 
time. 

The second model assumed two dimensional heat transfer 
(radial and circumferential). The same physical model was 
used as in the one.dimensional case with the sole exception 
that 48 nodes were used instead of 12. The representative 
nodal points and an example of the thermal connections from 
one of the nodal points are shown in Figure 28 in Appendix C. 
The four nodes on the surface of the storage container were 


modeled as zero volume boundary nodes. The sinusoidal and 


tabular representations were used to describe the surface 
temperature of the storage container at each boundary node. 
Actual data taken at each position, rather than bulk data, 


were used as the input data for these representations. 
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The same assumptions made in the one dimensional case 
were also applicable to the two dimensional model. A com- 
plete discussion of the calculation of the radiation ex- 
change factors in the two dimensional case is given in 
Appendix C. 

The effect of natural convection was studied in both 
the one and two dimensional models. References 5 and 7 
give correlations between the Grashof number based on the 
gap width and the effective thermal conductivity. The 
Grashof number was calculated from the equation 


p*gB(AT) 6° 
G SS eg 
x 2 
where 6 = width of the air gap 


AT = maximum temperature difference at any instant 
of time in the air gap 


B = : where T = 565°R 


At T = 565°R, air has the following properties 

p = 0.07 1bm/£t? 

yw = 0.046 Ibm/hr-ft 
The maximum Grashof number for this experiment was calcula- 
ted to be 1.25 x 10°, The diameter ratio was approximately 
1.5 and the log Gr = 6.1. From the Beckmann correlation 


k 
[Ref. 5], the effective thermal conductivity ratio ( was 


approximately 3.2. Using the Liu correlation [Ref. 7] 


k. Pr” Gr 5 0.278 
bE 7 9-435 | ingetpr oe 


where the Prandtl Number = 0.707. An effective thermal con- 


ductivity of 4.0 was assumed as the average value of these 
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two correlations and it was used to study the effects of 


free convection. This change was placed into the TRUMP pro- 


gram by increasing the value of the thermal conductivity of 


‘ air by a factor of 4 in each of the TRUMP runs. 
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V. RESULTS 
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A. ANALYTICAL MODEL 
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. Using the sinusoidal temperature distribution as an ap- 
proximation to the actual average experimental data as 
shown in Figure 7, comparisons were made between predicted 
temperatures and actual temperatures for two radial loca~ 
tions in the rocket motor. Figure 12 compares the results 


on the surface of the rocket motor and Figure 13 does the 
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comparison at the center of the rocket motor. An uncertain-~ 


Tren ie 


' ty analysis is given in Appendix E which establishes the 
uncertainty bounds for both the predicted and the actual 
temperatures. These uncertainty bounds are included in 


; Figures 12 and 13. 
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It is readily seen from Figure 7 that a sine wave was 
not an ideal fit as an approximation to the experimental 
data, as it varies as much as 20°F during part of the cycle. 
However, it was also noted that the sine wave closely approxi- 
mated the experimental data during the heating phase of the 
! cycle and: only during the cooling phase were there large 
{ variations. As the main purpose of this study was to design 
a model that would be useful in optimizing storage container 


a design, the errors in the cooling phase are not critical as 


long as the temperatures reach the same minimum point before 
r beginning another cycle. Figure 12 shows that the maximum 
surface temperature of the rocket motor predicted by the 


analytical model is a good approximation to the actual 
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experimental data. Again it is noted that, ‘in actuality, 
the motor cools faster than the predicted value. The maxi- 
mum difference in temperatures on the surface of the motor 


. fs 15°F, Figure 13 shows that the predicted value and the 


experimental value of the temperature at the center of the 
rocket motor were in close approximation except during the 
early stages of the cooling phase where a maximum tempera- 


ture variation of about 5°F occurred. 
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One of the reasons the system cools faster than predicted 
could be the Jight breeze that is usually evident in the 


early afternoon hours at China Lake that is not present dur- 


from the wind to study the effects of a light breeze on the 


nsec Nie cal ma a 


surface temperature of the storage container. 


i 

ing the morning. No attempt was made to shield the system | 
Another point not taken into account by the analytical 
model is the fact that the time delay at any point in the 

system is not constant throughout the day as predicted in | 
Figure 10, but varies as given by the data in Appendix D. 
Time delays between the peak temperature on the container 

wie Pace and the peak temperature at the center of the rocket | 

motor vary from about 250 to 400 minutes, whereas the low I 

i 

temperature on the surface of the container and the low ten- : 

perature at the center of the rocket motor vary from about 

150 to 250 minutes. The analytical model predicts a constant 


. variation of 388 minutes at the center of the rocket motor 


and 159 minutes at -tthe surface. 


Fae te 


B. TRUMP MODEL 
1. One Dimensional 
Four variations of the one dimensional TRUMP model 


. were investigated and compared to the experimental data. 
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Figures 14 and 15 compare the TRUMP predictions to the act- 
ual experimental data at the surface and the center of the 
rocket motor,respectively. The TRUMP variation used for 
this comparison modeled the storage container temperature 
with tabular ents (See Figure 11) aud-aseumed convection 


was present (— = 4.0). The uncertainty analysis (Appendix 
k 


E) established the uncertainty bounds for both the experi- 


| 
| 


mental and the analytical data in these Figures. The varia- 
tion between the bulk temperature predicted’ by TRUMP and 

the experimental data ciosely matches with only two experi- 
mental points in Figure 14 falling outside the uncertainty 
bounds for this one dimensional model. Figure ll shows 

that the tabular data that TRUMP interpolates is a good 
approximatton to the averaged experimental data. At ‘the 
center of the motor, as shown in Figure 15, all experimental 
points fall within the predicted error bounds. A comparison 
of the four one-dimensional TRUMP variations are given in 
Figures 16 and 17 at the surface and the center of the rocket 
motor respectively. It is clearly seen from these Figures. 
that the convection assumption results in.an increase of 2°F 
in the maximum temperature and a decrease of 2°F in the 
minimum temperature on the surface of the rocket motor. This 


temperature change drops to + 1.5°F at the center of the 
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motor as shown in Figure 17. The differences between the 
sinusoidal approximation anc the tabular approximation of 
the experimental data was clearly shown in Figure 11. The 
data in Figures 16 and 17 can be easily correlated to that 
in Figure 11, thereby explaining the differences in the pre- 
dicted values. 
2. Two Dimensional 

Four variations of the two dimensional TRUMP model 
were investigated and compared to the experimental data. 
Comparisons of each TRUMP variation to the experimental data 
are given in Figures 18 and 19 for node 8 (located on the 
skin of the rocket motor at the 1200 position) and node 1 
(at the center of the rocket motor) respectively. These 
Figures show that the TRUMP variations that used tabular 
data to model the surface temperature of the storage con- 
tainecr predicted temperatures that more closely approximated 
the experimental values than were those predicted by TRUMP 
variations using sinusoidal data to model the surface ten- 
perature. Appendix D shows that all points on. the surface 
of the storage container reach their minimum temperature at 
the same time; however, these points reach their maximum 
temperature as much as 200 minutes apart. Whereas, all the 
points on the surface of the storage container are in -:phase 
at the minimum temperature, they rapidly become out of phase 
as the container temperature rises. This varying phase 
shift makes it difficult to model the four boundary nodes 


with sinusoidal approximations which must have constant 
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phase shifts. Sizable errors in the input data during some . 
{ 
q 
; | 
Zz 
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parts of the cycle were caused by these varying phase shifts. 

These errors in the input data led to the variations in the 

predicted temperature values. As noted in the one dimen- 

sional section, the inclusion of convection effects does 

not produce large variations in the predicted temperatures. 
Figures 20 and 21 show the actual temperature dis-~- 

tributions on the surface of the storage container and on 

the surface of the rocket motor respectively at maximum 

bulk temperature compared to a two dimensional TRUMP program. 


The TRUMP variation used for this comparison assumed no free 


convection in the air gap and used tabular data to approxi- 


! 

Mate the surface temperature .of the storage container. 
C. GENERAL 

A comparison was made between surface temperatures on 
the storage container that contained the rocket motor and 
the storage container that was empty. “The low temperature 
was about 4°F colder in the empty container, whereas the 
high temperature was about 4°F higher on the container that 
contained the rocket motor. The empty container had a 
faster response time than the one containing the motor. i 
The differences in heat capacities, radiation effects, and i 
natural convection all contribute to the changes in tempera- : 


ture noted. 


Figure 20: Temperature Distribution at Surface of Storage 
Container at Maximum Bulk Temperature at approximately 
1500 on 2 August 1972. 


FEDER WEEN COT OES LANE OSE WENT DIO ETON EE eI? 


meets cee natal A crete se, 


1200 


~-- Actual 
0 2-D TRUMP 


0600 


Figure 21: Temperature Distribution at Surface of- the 
Rocket Motor at Maximum Bulk Temperature -at approxi- 
mately 1500 on 2 August 1972. 
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D. LIQUID CRYSTALS 


Shearer neem 


the surface of the storage container gave brilliant colors 
- under the intense desert sun. These colors were much 


clearer and brighter than the same crystals viewed under 
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laboratory lighting conditions. The liquid crystals phvto- 
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The encapsulated cholesteric liquid crystals applied to 
graphed well in both the color movies and the color slides. 
The movies showed by time lapse photography the rapidly 
changing surface temperature of the storage container. Two 


sample color prints maie from the color slides are enclosed 


t 
[ 
| 
| 


as Figures 22 and 23 to show the brilliance of the colors 
and the feasibility of obtaining data from color photos. 
The only photographic problem encountered was the intense 
reflection of the sunlight from the polyurethane film. This 
problem was partially overcome by taking the photographs 
from angles where the reflection was less intense. Quali- 
tatively the liquid crystals were not adversely affected by 
the sun's rays after two weeks of desert exposure. No 
accurate quantitative determination was attempted; however, 


rough approximations were made at the site. These approxi- 


mations were made by noting the color exhibited by a crystal 


at a certain time and then comparing the calibration. of the 


crystal (Table 1) to the temperature recorded by the thermo- 


| 
: 
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couple located directly beneath the region of color change. 
Z The readings were within + 2°F, which was very encouraging, 
especially considering the approximations made while taking 


these measurements. Although photos were taken only during 
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the initial two weeks of the study, on site observations in- 


dicate that the cxystals are still showing brilliant colors 
after 3-1/2 months. Preliminary onaenie indicates that the 
polyurethane film did protect the crystals from decomposing 
from the sun's rays and from being worn away by the wind 
blown sand. 

It was noted that the surface temperature of the storage 
container under the liquid cyrstals reached temperatures up 
to 15°F higher than a similar point not under the crystals. 
This 15°F difference was only evident when maximum tempera- 
tures were obtained, During sunlight hours ‘the temperature 
under the liquid crystals was always somewhat higher; how- 
ever, at night both temperatures were equal. The difference 
in the container surface temperatures led to a difference of 
4°F on the surface of the rocket motor and 1°F at the center 
of the motor. it is believed that the difference in emissi- i 
vities of the gray and black surfaces resulted in the differ- 


ence in container surface temperature. 


VI. CONCLUSIONS 


From the results of this investigation, the following 
conclusions were drawn: 

1. Although a sine wave is not a perfect fit for the 
experimental data at all points, it is useful in predicting 
bulk temperatures in the rocket motor, especially if only 
the high and low bulk temperatures are of concern. This is 
especially true in the one dimensional case. In the two 
dimensiona]. case, the problem of phase shift variations 


make the method of sinusoidai variation less desirable, 


~~ 


although still useful. 

2. The simulation of the actual data by a table of 
. temperatures gave the most accurate predictions of the 
experimental data. This method should be used whenever tubu- 
lated data are available; however, this gilt penevaliy not 
‘be the case for design work, in which case the sinusoidal 
approximations must be used. 

3. The flexibility of both the analytical and computer 
simulations allow the changing of many parameters. The re- 
sulting effects of these changes on rocket motor temperature 
can ‘be studied with the models. 

4. The convection assumption for this system resulted 


in only small changes in temperature and can be neglected 
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when predicting design temperatures. Either the Liu or 
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Beckmann correlation should be used to determine if convec~ i 
/ 
tion can be neglected in a particular system. | 
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5. The use of an empty storage container to obtain sur- 
face temperature data is a good approximation to vsing one 
with a rocket motor inside. 

6. It is feasible to use liquid crystals for thermal 
mapping under desert conditions. Color photography with 
standard equipment gives excellent results since brilliant 
colors were observed. 

7. ‘The liquid crystals appear to be stable for at least 
two weeks under the desert conditions when protected with a 
polyurethane coating. 

8. The application of the liquid crystal system to the 
surface of the storage container resulted in-large increases 
in the surface temperature of the container throughout the 
hottest part of the day. Care must be taken in applying 
and interpreting thermal readings from liquid crystals when 


exposed to radiant heating. 
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VII. RECOMMENDATIONS 


From the results of this basic study, the following 
ecoohmendations for £uture work are offered: 

1. To refine the results of this project, a second 
rocket motor storage container system should be instrumented 
with the following changes: 

a. Liquid crystals should not be applied to the 
system used as the experimental neaed: As steel is a good 
thermal conductor, axial conduction on the surface of the 
storage container may be significant. Heat flow from the, 
area where the crystals are applied may lead to higher than 


normal temperatures at other points on the surface of the 


container, 


after the loading of the dry sand so that‘an accurate deter- 


mination of the density of the propellant simulant can be 


Lees 


! 
| 
| 
| b. The rocket motor should be weighed before and 
| 
{ 
| determined, 

| c. Four additional thermocouples should be located 
| on the surface of the storage container to better enable 

the averaging of data. At present, the #1 thermocouple 
which was used as the average temperature reading of the 

top quarter of the surface of the container, in actuality is 
its hottest point; likewise the #10 thermocouple was used 
as the average temperature of the bottom quarter of the 


surface, in actuality it's the coldest point. For averaging 


data, it is recommended that thermocouples be placed at 0130, 


| 
i 
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0430, 0730, and 1030 and the quarters of the system be di- 


vided at 0300, 0600, 0900, and 1200 to give a more realistic 
bulk temperature. Thermocouples at 1200 and 0600 will pro- 
vide the maximum and minimum temperature of the system. 

2. The TRUMP program should be rerun in both the one 
and two dimensional form, varying the mesh sizes to deter- 
mine the optimum number of nodes. 

3. A long term study of the effects of the desert en- 
vironment on liquid crystals should be done. The crystals 
should be calibrated before being placed in the desert and 
then brought to a laboratory for recalibration at specific 
intervals. 

4. Several modifications should be made to the TRUMP 
program to make it comparable to the version used at Law- 
rence Radiation Laboratory. The variable conductivity sec- 
tion (BLOCK 2) aud the PLOT subroutine (BLOCK 11) need to 
be corrected. The TIMEP subroutine which allows the setting 
of the problem time interval between data output should be 
added to this version of TRUMP. It would also be advanta- 
ae increase the amount of tabular data that could be 
read in as boundary temperatures, 

5. Erom an academic standpoint, the effects of free 
convection in an air gap with varying boundary temperatures 


should be investigated. 
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APPENDIX A 


Introduction to Liquid Crystals 


Liquid crystals were first discovered in 1889 by Reinit- 
zer [Ref. 8] and the favesminekiens of Lehmann which con- 
tinued to 1915. Liquid crystals were considered to be 
laboratory curiosities with no scientific or practical merit 
until the 1950's. They share some of the properties of both 
liquids and eeyaeeies for example, a typical liquid crystal 
Substance scatters light in symmetrical patterns and re- 
flects different colors depending on the angle from which 
it is viewed. Studies in the last few years have helped to 
clarify the unusual molecular structure of liquid crystals. 
Many applications arise from their ability to detect minute 
fluctuations in cemperature, mechanical stress, electromag- 
netic radiation and chemical environment by changes in their 
color. 

Liquid crystals are divided into three classes; smectic, 
nematic, and cholesteric, depending on the degree of spatial 
arrangement of the molecules in the mass of the material and 
the type of the material [Ref. 9]. In this project only 
cholesteric liquid crystals were used and therefore only 
their properties will be mentioned. The molecular structure 
of cholesteric liquid crystals is characteristic of the 
esters of cholesterol (Figure 24). The molecular layers 
are very thin with the long axis of the molecules parallel 


to the plane of the layers. The individual molecules are 
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Figure 24: Molecular Structure of Cholesteric Ester 


: | 
; ; ; AY viewes 
incident bean retlected beam F 
liquid crystal fiin af transmitted beam : 
2 , i 
« 1 ‘ 

' 


TA eee 


ed men! see wt any 


ff 


black background absorbed bean- 


Figure 25. Light Reflection from Liquid Crystals. 
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basically flat, with. a side chain of methyl groups (-CH,) 
projecting upward from the plane of each molecule. This 


configuration causes the direction of the long axis of the 


molecules in each layer to be displaced slightly from the 
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corresponding direction in adjacent layers, This displace- 


ment, which averages about fifteen minutes of are per layer, 
ts cumulative through successive layers, and the overall 
displacement traces out a helical path. 

The molecular structure of cholesteric liquid crystals 


gives rise to many peculiar optical properties. If linearly 


~ hws 


‘polarized light is transmitted perpendicularly to the molec- 


Tight will be rotated to the left in a helical path. There- 
fore, the plane of polarization will also be rotated to the 
left, through an angle proportional to the thickuess of the 


transmitting material. Liquid crystals are the most opti- 


ular layers, the direction of the electric vector of the 
t 


cally active substances known, Another strictly crystalline 
sokdeaa property exhibited by cholesteric liquid crystals 
ts circular sidice ten: When ordinary white light is inci- 
dent to a cholesteric material, the light is separated into. 
i two components, ene with the electric vector rotating clock- 
wise, the other rotating counterclockwise. Depending on 
the material, one of these components is transmitted, and 
; the other is reflected. It is this property that gives the 
| cholesteric phase its iridescent color when it is illumina- 
ted by white light. The particular combination of colors 


depends on the material, the temperature, and the angle of 


the incident light. 


The molecular structure of a cholesteric substance is 


very delicately balanced and is easily upset. Any small 
disturbance that interferes with the weak forces between the 
molecules can produce marked changes in optical properties 
such as reflection, transmission, birefringence, circular 
dichroism, optical activity and color. The must striking 
optical transformation that occurs in a cholesteric sub- 
Stance, in response to small changes in its environment, is 
the variation of color with temperature. The crystal lat- 
tice is disrupted by the thermal vibrations giving succes- 
sive transitions between the solid, the mesophase, and the 
isotropic liquid with rising temperature. The change from 
the three dimensional order of the crystal lattice to the 
disorder of the isotrepic Jiquid occurs via one or more in- 
termediate states, each of which has a particular tempera- 
ture range at which it is stable [Ref. 10]. | 

A cholesteric liquid crystal system responds to changes. 
in temperature by sequentially passing through the complete 
visual spectrum (red through violet) in fractions or multi- 
degrees, depending on which cholesterol esters comprise the 
formulation. This color phenomenon is reversible and has 
been reported to function over a temperature range of ~20°C 
to 250°C. A very important point to note is that at a cer- 
tain temperature a given material or combination of materials 
Will always exkibit the same solor. Also, the rate of change 
from color to color as well as the exact temperature at 


which the specific color changes -occur are invariable. By 
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desired temperature combination can be obtained. The thick- 
ness of the cholesteric film does not affect the predominant 
wave length of the reflected light; the light becomes cir~ 
cularly polarized [Ref. 11]. 


The colors scattered by the liquid crystals represent 


| 
| 
: 


only a fraction of the incident light (Figure 25). The re- 
maining portion of the incident light is transmitted by the 


liquid crystals. Therefore, an absorptive black background 


| 
mixing cholesteric substances in various proportions, any 
| 


must be used to prevent reflection of the transmitted light, 


thereby enhancing the resolution of the scattere’ colors or 
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wavelengths reflected by the liquid crystal system. 

The cholesteric liquid crystal systems often present a 

r number of problems due to the fact that they are viscous 

liquids. Some problems associated with the handling and 
the use of these materials are: 

1, The tendency of the liquid crystal system to flow 
during application can cause vautatione in applied film 
thickness. This may result in non-uniform thermal patterns. 


2. Direct exposure of liquid crystals to adverse en- 


a nace ee le Ee Vee 


vironmental effects can cause variations in their sensiti- 


a bet 


a 


vity and deteriorate their color response in a few days. 
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| . These problems can be partially overcome by using an 
\ encapsulated liquid crystal material system. The capsules 
. he. ‘are 20-30 microns in diameter and are a water-based slurry 


suitable for application by conventional coating techniques 


such as brushing or spraying. 


Encapsulated liquid crystals offer several advantages: 
1. They convert the liquid crystal system to a pseudo- 
solid, which provides for easier handling, application, and 
; use, 
2. They provide longer shelf life by minimizing surface 
contamination and giving protection from ultraviolet light 


[Ref. 12]. 


3. They exhibit relatively unlimited fatigue life. 


4. They reduce the angular dependence of the color 


observed. 
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APPENDIX B 


Analytical Solution 


The method of complex temperature as presented by Arpact 
[Ref. 1] was used to find the steady periodic solution of a 
body experiencing a periodic sinusoidal disturbance. The 
general heat conduction equation in cylindrical coordinates 
was the basis for this derivation. It was assumed that no 
heat sources existed in this problem, that the rocket motor 


storage container system was infinitely long, that there was 


no heat conduction in the axial or circumferential -direc- 


tions, and that the container surface temperature was spa- 


cially uniform. Figure 26 gives a basic sketch of the system. 


The assumptions reduced the heat conduction equation to 


oT 
dae Pe = 1 of (X) 
Yr or Qa ot - 


where r is the radial distance from the center of the rocket 


motor, T is the temperature of the rocket motor at time t 


and position r, and a is the thermal diffusivity, a property 


of the conducting material. 


Paw : 
GF as (2) 


where k is the thermal conductivity of the conducting 


material, p is the density of the material, end ec is the 


specific heat. All thermal properties were assumed: to be 


eres v- ie we er 


constant over the temperature range of this problen. 
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Figure 26. Analytical Model of Experimental Systen. 


The boundary conditions used in this derivation were 


at = 0 at r= 0 
QT h = 
and aor te (I - T,) at r= to 


where ro is the inner radius of the rocket motor. 
T,, is the known storage container temperature which is 
assumed to vary as 


T= (Iy- T,) sin wt + T, 


Ty = maximum bulk temperature of the storage container 


A ™ average bulk temperature of the storage container 
27 
) 


T 
Ww = frequency of the sinusoidal variztion GEhoars. 
: 2 


= time 
h is the effective heat transfer coefficient across the 
air gap between the storage container and the rocket motor. 
It combines the heat transfer effects of radiation, convec- 
tion, and conduction into one coefficient, The radiation 
-coefficient was linearized by assuming constant temperatures 
(Ty > To)» representative of the average temperatures expec~ 
ted in the problem, in the equation 
bray = 71-297, #1) (7 + 75°) 

where ¢ is une Stefan-Boltzmann constant and # 


1-2. 


radiation exchange factor between surfaces 1 and 2. The 


is the 


convection coefficient is 


where k. is the effective conductivity of air as obtained 
from the Beckmann and Liu correlations [Ref. 5 and 7] and 


6 is the width of the air gap. In the analytical model, 
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the effective conductivity was assumed to equal the conducti- 


vity, thereby treating it as pure conduction and 


Equation (1) was non-dimensionalized using the following 


relationships 


T-T 
6: = (a non-dimensional temperature) 
T,- Th 


E = = (a non-dimensional distance) 
o 
to give 
00 2 
1 dCEF) = Fo ae 
E gE O t 
with boundary conditions 
de = 0 at €=0 
. AO 
and JF =-8(0 - sin wt)at § = 1 
hr | 
where ff = 5 Sl is the Biot modulus. (which compares the 


relative magnitudes of the effective heat transfer -coeffi- 
cient across the air gap and the internal conduction 
‘yesistances to heat transfer): 5 
An initial condition was not specified as the only 
concern was with the steady state, periodic behavior. Fol- 
lowing Arpaci [Ref. 1], a complex temperature was defined as 
W(r,t) = 6° (r,t) + i6(r,t) 


where (r,t) satisfied 


aceddy xr . a | 
Ee a ae @) 
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a ea 
| with boundary conditions 
: : BoP +1 Bo at E = 0 | 
: and wy = B62) e -8(6*-cos wt)+ 1{-BeO~sin wt)}at E = 4 
This leads to 6(r,t) which satisfied | 
H 
| 2 ep ; Fo do 
| with boundary conditions 

and 4 = -B(6 - sin wt) at Es] 


also 6" (r,t) which Satisfied 


with boundary conditions 
* 
fy = 0 at € = 0 
* 
“ = ~B(0*~cos Wt) at Ex 


and 


A solution of the form 


te ae a ee een 


W(ryt) = $(x) T(t) 


was assumed, where for large values of time t(t) was 


assumed to. equal etut, therefore, 


W(r,t) = o(r) ett (4) 
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Equation (4) was then substituted into equation (3) 
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with boundary conditions 
“et = 0 at & 


and ae = ~B(o-1) at Fel 


Wt 
oO 


Equation (5) was expanded to give 


2 | 


2 iwr 
ale 4,2 go. 2. = 0 (6) 


t& 


Now, let Z = 


and substitute into equation (6) 
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with boundary conditions 


ae =0 | at Z = 0 


fier, 
—— (6-1) at Z = a 
°o 


The general solution of equation (7) is 


roy 
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and 


»n 
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as given in Ref. 13 with 


i) 

3 
ete here eee Stl aeons ARR RE te et tm nn 
peeks Meaticat a riled! te 


1,.4 oo 4 Bene 
(52) . (52) i 
1,(2) = 14(32)7+ 2 + ---- = ) ees : 
‘ | (21)* (nt) 
n=0 i 
and . bz) 20 
Ky(2) = ~{ytlo(32)}1,(2)+ DY at p+ Gt B 
~~ 


Now using the boundary condition 


ao 0 at Z = 0 
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and differentiating equation (8) yields 


3 d(1_(2)) a(K (Z)) 
st = ¢ ° +c 2 
: z 1” dz 2° dz 
where 
Py (I (Z)) 
coer \ eames at Z = 0 
and 
d(K, (Z)) 
oo # 0 . at Z = 0 
therefore Cy = 0 
and ¢= C, 15(2) (9) 


Now using the second boundary condition that 


{ 
{ 
4 
; 
4 
j 
| 
| 
: 
! 


Lwr 
d 
Soe - —B— (g-ay at 2 gf (10) 
fiwr 
° 
a 
and differentiating equation (9) gives 
A(z, (2)) 
dg = ¢ 9 
dZ 1 az 
(I, (2)) 
Noting that age I, (2) and substituting into 
equation (10). 
B ( iwr * 
Cc C,l \@——— /-1 
a 


Rearranging and solving for Cc, 


then substituting into equation (9) 


(2) 
6 = (11) 
( iwr ry ieio , Wz a 
I t= 
Wat) +3 ) 
Now as 
J, Cimx) = 1°T, (mx) [Ref£. 3, p. 135] 
* Ey 2 B, 2 
Wr wr 
q ° 1/25) * ( o 3/2 ) 
I, ( a i tj = Je aa i gE 
and 
wr 2 ( wr 
o ,1/2 ,)_ 2 o .3/2 ,) 
Ty ( ria a ale a | co RY 
wr 2 
Let a = a and substitute into equation (11) 
J (3/2 se) 
a ea a (12) 
J_.(4°'°" ajo Si J rico a) 
re) B 
3/2 _ 49% 31 3m 1 
Now i = et “4 = cos ok + i sin “ S (-1 + i) 
v2. 


Substitut“ 1g this into equation (12) 


> nd ad) (13) 
3/2 3/2... 
J_(i y+ — (1- -i)3, (i a): 
6 a 5a i 


Say7 ish, 
As J (abi )= J, (abe 4) = BER (a&)+ iBEi (aé) 


and 311 
3, ¢ab43/*)= J, (abe™ 4) = BER, (a&)+ iBEi, (af) 


Substituting these results into equation (13) yields 


BER, (a&) + iBEi 9 fab) 


BER (a)+ iBEi 9(a)+ a5 (1- 4) (BER, (a)+ iBEi, (a)) 


(14) 


eee ee 
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After rearrangement: 


BER, (aé) + iBEi (aé) 


b rr reece ee EN A y ooneeneeeaeennnemeneleneenenanetl 
[BER (a)+ —*-nER, (a)+ BE, (a)}+ 4 [BEi, (a)+ BEI, (a)- —2-BER, (a)J 
Letting v2p V2B 2 ¥2B 
X. = BER (a)+ —*-BER,(a)+ —*-BEi, (a) 
R fe) 2B. 1 /28 1 
and 


a 


xX, = BEL (a)+ 


_ da : 
, BEi, (a)~- Fe 


5 
% 


. . .and substituting into,equation (14) gives |. 


BER (af) + iBEi (af) 


Xp + ix, 
Oe tee ww a heshonslizing..che denowinaton ie Jas.. ~ om es peu oe ek pe te nate 
BER (ab) + iBEi, (ab) 
= — — — ~ (X, - ixX,) (15) 
xX 2 +X 2 R 2 
R i 


Now 


ie BER, (a6) pABEL, (a6) X,) +4 (BEL (ab) Xp BER, (a6)X,) 
a + af 


which after rearrangement gives 


we oe BEi “(aby gg 


jh dina nds we estue hd Cara dentate tina area ra cont ae athe tae tai: 


og e (16) 
Xp + x, 
where 
, _, BEA (a&)X, - BER (a&)X. 
Re bag 1 ° R Ce) i 


BER, (ag) X, + BE, (ab)xX, 
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BER “(a&)+BEi_ “ (ag). ¥ e 
Ww(r,t) = So [cos(wttd )+i sin(wtt+6é )] 
: x 


See ee 


. Som wee BER: 5 Cab) ¥BEL, Beary: ee ee ee i ee TT) tele nem ee ee 
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Substituting into equation (4) gives 


oer 2 ea on 
{BER (a&)+ BEi, (a&) 
2 2 


Xp: + Xy 


* 
ei (wt +6 ) 


Wr, t) 


which also equals 


Rt % 
As this problem was modeled as a sine wave, the imaginary 


part. of (r,t) was nee 


which is the analytical solution of infinitely long -concen- 
te et Ae Be ee te tee ee ce ewe me He eee Fhe web te 8 ete eet sete ew ae fee 
tric cylinders experiencing a periodic sinusoidal temperature 
variation on its outermost surface when “eat conduction is 


assumed to be radial only-.. 


In summary 


| 
T-T [BER (aE) +BEi~(aé) , 4 
‘ A Cee o . * 
O(r,t)=— =F Sg ean. Cte) 
MTA X,” + X; j 
R i | 
where | 
T = the temperature of a point r in the rocket motor at | 
time t i } 
T= average bulk temperature of the storage container f i 
he Maximum bulk temperature of the storage container , i 
= frequency of the sinusoidal variation(27/24 hours) : F 
t = time ‘ 
§ j 
Es = = dimensionless distance from the center of the rg 


Xo rocket -motor 
xr = distance to the surface of the rocket motor 


° 
x = distance from the center of the rocket motor : 

fur, - 
a= a conduction parameter ; 
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Qa = De 


= thermal diffusivity 


density 

thermal conductivity 
specific heat 

real Bessel Function 
imaginary Bessel Function 
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X, = BER (a) + —+ BER. (a) + —* BEt. (a) 
re) /28 1 /78 1 


X, = BEi (a) + —\ BEi,(a) - —* BER, (a) 
fr ° /ze * jig =? 

B = ~~~ = Biot modulus 

st 7 ean72 BEi (ab) Xx, - BER, (a&)X, 
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The following computer programs were used to investigate 


a wide variety of parameters, The outputs are samples of 
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APPENDIX C 


TRUMP Solution 


TRUMP is a computer program for solving transient and 
steady-state temperature distributions in multidimensional 
systems. This program was developed in 1965 at the Lawrence 
Radiation Laboratory by A. L. Edwards [Ref. 2] for their 
cDC/3600 computer. The program was adapted to the Naval 
Postgraduate School IBM/360 Model 67 computer system in 
1971 by C. Erbayrum [Ref. 3] from a version used by the 


B. F. Goodrich Corporation. 
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TRUMP is a mult’ durpose program able to solve a wide 
variety of problems involving flow in various kinds of 


potential fields such as heat flow in a temperature field. 


TRUMP allows the solution of general nonlinear parabolic 
partial differential equations both in steady-state and 
transient problems. Complex geometric configurations with 
multidimensional flow may be solved using various coordinate 
systems. Initial conditions may vary with spatial position. 
Material properties, boundary conditions, and other problem 
parameters may vary with spatial position, time, or the 
primary dependent variable, ° 


Input data are fed to TRUMP in "Block" form through its 


12 input data blocks, A complete description of each of 


* 


these blocks is gi en in Ref. 2. A model of the problem 


must be constructed and data from this model read into 


TRUMP through the data blocks. 
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Two models were used to simulate the rocket motor stor- 
age container systz2m and several variations of each model 
were investigated. 

The first model assumed one dimensional heat transfer 
(radial) with the assumptions that the system was infinitely 


long and that the container surface temperature was spatial- 
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ly uniform. The system was modeled as two infinitely long 
concentric cylinders separated by a 2.94 inch air gap. The 


inner cylinder was constructed of 4130 steel and was filled 


materials used in the experimental system are given in 


SFR ET PT LPI Pe TT ETI SITY TE LR PT TT 
‘ 1 ; 4 i 


| 
{ 
i] 
| 
| with dry wind blown sand. The thermal properties of the 


Table IL with units most easily compared to the actual data 


obtained from the system at China Lake. 


TABLE IL 
Thermal Properties of Materials 


Material Density 3 Specific Heat Thermal Conductivity 
Sand 0.05486lbm/in~ O.195BTU/1bm°F 0.00026BTU/min-in-°F 
Steel 0.2807 1bm/in3 0.109BTU/1lbm°F 0.364 BYU/min-in-°F 
Air ‘0.0000436 " 0.240BTU/1bm°F 0.:0000225 - 


with the representative nodal points given in Figure 27. 


mre a ae Wa oe 


The model was subdivided into volume elements or nodes 
Although the representative nodal point may be located any- 

where in the node or on the surface of the node, in transient 

problems it is usualiy located so that the lines connecting 

the nodal points are perpendicularly bisected by the con- 

nected area. This gives maximum accuracy. Two boundary 


conditions were given to the surface node. The first was a 


sinusoidal disturbance which closely modeled the actual 
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Figure 27. 
Model, 


‘Location of Nodes for One Dimensional TRUMP 
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average experimental data obtained at China Lake as‘given in 
Appendix D. The sine wave exhibited a maximum temperature 
of 138°F and an average temperature of 104°F. Its period 
was 24 hours (1440 minutes). The second boundary condition 
was the actual average surfuce temperature of the storage 
container given at two hour intervals. Both these boundary 
conditions are approximations of the actual] surface tempera- 
ture. Two hour intervals were the minimum allowable for the 
tabular data as this version of TRUMP has a maximum table 
size of 12. 

Several assumptions were initially made. The thermo- 
couple data obtained from the experiment at China Lake gave 
the average temperature at a point on the storage container 
and not the actual outside surface temperature. As this 
container wall was only 1/16 of an inch thick and made of a 


good thermal conductor, it was decided to model this data as 


a zero volume boundary node with a known temperature impressed 


on it. It was also assumed that heat transfer across the 


air gap occurred only by radiation and conduction, neglecting 


the effects of free convection. 

It was estimated that the surface emissivities for the 
rocket motor and the storage container were 0.9 [Ref. 6]: 
based on their haze gray surfaces. The radiation exchange 


factor for this geometry is given by 
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A sample input deck for the tabular approximation of the 
boundary condition is given at the end of this appendix. 
Several cycles of output data for the one dimensional model 
are also given, 
The second model assumed two dimensional heat transfer 
(radial and circumferential) with the assumption that the 
system was infinitely long. The same physical model was 
assumed for the system except 48 nodes were used instead of 
12. The representative nodal points are given in Figure 28. 
The four surface nodes (12, 24, 36, and 48) each had two 
F different temperature approximations applied, a sinusoidal 
representation and a tabular input taken at two hour inter- 
vals. The four surface nodes were also modeled as zero 
volume boundary nodes. Each internal thermal connection 
between nodes is described in the input data by specifying 
the two node identification numbers, two connector lengths, 
and two interface dimensional factors. An example of the 
thermal connections of node 4 is shown in Figure 28 and the 
input data in BLOCK 5 of the two dimensional TRUMP program. 
The calculation of the radiosities in the two dimension-~ 
al case Wau aecoupiicned by using a radiation-network and 
the method of crossed-strings. 
The radiation shape factors for the two dimensional sys- 


tem were determined by the method of crossed-strings [Ref. 


14}. The graphical construction for this method is given 


in Figure 29, 
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Location of Nodes for Two Dimensional TRUMP 
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Figure 29: Graphical Construction for Crossed-Strings 
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Eon is defined as the fraction of energy leaving sur- 


face m which directly reaches surface n. From the physical 
dimensions of the model A,= A,= Ac= 9.45 59:42 ye alm A 


in., 3 5 
. = 14.05 a Ant = 12.65 Sq-2n assuming unit depth. Let 


Sy equal the length of A 


; 
: 
I 
7 
E 
t 
: 


i° 


The crossed-string method lets each surface represent 


the effective surface obtained by stretching a string tight- 
ly over the radiating face between the bounding edges, to 
produce a surface that cannot see any of itself. For ex- 
ample, surface 2* in Figure 29 stretched over surface 2. 


By definition Foy = 1, which by reciprocity leads to 


Ay 
3 Foor = oa Fora = 0.9 
; and therefore since 


To calculate the direct'radiant heat exchanged between 
surfaces 1 and 2, a minimum-length line was stretched con- 
necting edge B: of A, to edge E of A, and a second minimum 
length line from edjze L of Ay to edge F of Aye These lines 


are labeled L, in Figure 29 and are equal to the width of 


1 


the air gap, L, = 2.9375 in. Minimum length lines were also 


Lt 
stretched from point B on Ay to F on Ay and L on Ay to E on 


Ay The length of these lines is Dd, and is made up of two 


parts; Dis the tangential distance from F to surface A, 
and D,", the arc length from the point the tangent hits A, ‘ 

. - Wo. ee 2 2 2. 6 E91 
to B, From geometry D, vr, ro 6.62 


1 oe = &.40" 
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therefore 4 
= ' oo " : 
dD, Di + Dy 11.04 
. 2D, - 2L i 
ae 1 ‘ 
Now Fi» amen -cracmaaa 0.86 
+ 3 
j 
From reciprocity, Ay Faia = Ay Foy 
x ; 
Oo { 
Fi. = —F = 0.578 | 
21 ry 12 
Now Fi3 is calculated from 
i 
Fi» + 2F 13 = J] 
therefore F = 0.07: 


13 : 


From symmetry Fi, = Fi3 


and then by reciprocity 


Ps 
; = 2 = 
Foy Bien Fi 0.047 


x 
Now Fo3 was caluculated y stretching minimum length lines 
from F to E, from E te G, from F to G and from E to E; where 


the length of the line from F to E = from E to G * So's 


from F to G = 2Dy> and from E to E = 0 
2s,' - 2D : 
2 1 a 
Fo3 = any Tae =» .113 
: i : 
F,. was now calculated from a4 
250 i: 
yy. PE ng P Fag. © F924. * Eps? . 
; fs 
therefore Fos = 0.002 L. 
As F,, was much smaller than the other F__, it was not Z. 2 
25 m-n Ef z 
included in the radiation-network diagram in Figure 30 H : 
[Ref. 15]. oe 
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Radiation Network. 
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; To calculate F,_, the blackbody potential of area 1 is 


set to unity and all other blackbody potentials are set as 


zero. An energy balance was written at each node giving a 


SR Teper po Tl 


. set of six simultaneous equations as follows: 
i Node 1 
[| Ay(U-Sq) = AVFy9(Sy-Jp)+ 24, Fy 3 (54 -F3) 
ae 
et Node 2 
boo : pes = 
> A, (0 Jo) AF, ,(J-5 pt 2A,F F545 ~I,)+ 2A5F 23649 -J 3) 
i 
f Node 3 ; 
: EA, (0-J5) = Ay F 12643 -JI,)+ AoFo3(5, ~Jo)+ AiF13053-J,)+ 
AiPy3(3g-Jg)t AlFy3 (3-55) 
Node 4 


A, (0-54) = AgFy,(Iy-Jo)+ AVFy4(J,-Jg)+ ADFo,(I4-Ts) 


Node 5 ? 
E,4,(0-J,) = 2A,F, 3(I5- J) t2AyFo,(I5-T,) +A, Fy 9 (I5-Jg) 


Node 6 
E,4, (0-Je) = 2A, F154 (Se-Ja)+ Ay Pia, ~J5) 


where Jn = radiosity of node n 


1 l-e Ay 


Eos = radiation shape factors previously calculated. 


Now the values of Ay» Ay and Fes were substituted into the 
energy balance equations which were then put into matrix 


form as shown in Table III. 
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Letting €,= Ey = 0.9, a standard computer solution for matrix 


problems gave the radiosities as listed in Table IV. 


‘TABLE IV 


Radiosities at Nodes 


Jy = .9046 
Jo = ,0529 
J, = .00495 
Jy = .000797 
J. = .000125 
J¢ = ,000080 


Now to find the radiation exchange factors from node 1 to 
nodes n, the radiation network shown in Figure 30 was re- 
duced to the equivalent network shown in Figure 31. Where 


the nodal equations are 


Ft 


1-2 A, (1-0) = EA, (J~-9) 
where E, hs 
= J = e7i 
1-2 A, 2 
134, (4-0) = E,A, (5370) 
E, A 
2 2 : 
F 3 * a 3 7 +066 


These values of the radiation exchange factor are used 
in the two-dimensional program. A sample input deck for 
the sinusoidal boundary condition is included at the end of 
this appendix. Several cycles of output data for the tabu- 


lar boundary condition are also given. 
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Equivalent Radiation Network 


Figure 31: 
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APPENDIX D 


Experimental Data 


The data presented in this appendix were obtained from 
the thermocouples on the rocket motor storage container sys- 
tem located at China Lake, California. The thermocouple 
output was read out on a Honeywell Electronik 25, 24 channel 
recorder which had been eatibeatad at 50, 100 and 150°F. 

The data was taken on two consecutive, typical summer days 
‘(August 1 and 2, 1972) at China Lake. Exch thermocouple 

was read once every 24 minutes. The first set. of data pre- 
sents the storage container temperature at four locations 
plus three different ways of averaging this data. It also 
presents the ambient temperature and the approximate time of 
day. The second set of data presents the surface temperature 
of the rocket motor and three ways to average this data. It 
also presents the temperature at the center of the rocket 
motor and the approximate time of day. Figure 32 shows the 
location of the thermocouples used. to collect this tempera~ 


ture data. 
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Figure 32: Thermocouple Locations for Experimental Data. 
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APPENDIX E 


Uncertainty Analysis 


An uncertainty analysis was carried out on both the ana- 


lytical solution and on a one dimensional TRUMP model of the 
rocket motor storage container system. In both models, the 
volumetric heat capacity of the sand (pC), the conductivity 
of the sand (k), and the emissivity of the surfaces were each 
varied by ten percent to determine the sensitivity of the 
system temperature response to each variation, Although 
other factors may also be varied, it was theorized that these 
three had the greatest effect on the heat transfer of the 
system, These factors were also known with the least accu- 
racy; the maximum uncertainty of each was estimated to be 
plus or minus ten percent (odds 20 to 1). 

In the analytical solution, varying the volumetric -heat 
capacity changed parameter a, varying the emissivity changed 
parameter §, and varying the conductivity changed both para- 
meters a and —. The effects on each parameter from each 


variation are given in Table V. 


TABLE V ? ee ag 
Change in Parameters due to Changes in Thermal Properties ; 
Change in Property Change in Parameters : 
Volumetric Heat Capacity + 10% at .12 : 
Volumetric Heat Capacity - 10% a-n- .12 
Emissivity + 102% 6B + .49 
Emissivity - 102 B - .37 
Conductivity + 102 a-~ .l1l, 8B - .22 : 
Conductivity - 102% a+ .13, B + .35 : 
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Each factor was varied holding the other factors constant. 
The changes in temperature and time delay were computed from 
the difference between these new values and those previously 
obtained From the analytical solut‘on. To obtain uncertainty 
bounds on the analytical curve, the second power equation 


[Ref. 16] was used, namely 


We = resulting uncertainty in the calculated temperature 
due to uncertainties in temperature caused by 


Wy = estimated uncertainty in volumetric heat capacity 


estimated uncertainty in conductivity 


& 
i) 


estimated uncertainty in emissivity 
An identical calculation was carried out to calculate the 
uncertainty in time delay. The results of these calculations 
are shown in Figures 12 and 13 for the Saveuee and center of 
the rocket motor respectively. The uncertainty in tempera- 
ture varied with time with a maximum variation of +2.75°F at 
the zenter of the motor and a maximum variation of +1.85°F at 
the surface of the rocket motor. The time delay varied by 
+31 minutes at the center of the motor and +11 minutes at the 
pated, The actual experimental data was also plotted on 
these Figures for comparison. 

The experimental data also had an uncertainty bound. 
Three primary factors made up this uncercainty bound; the 


accuracy of the thermocouple wire (i1.5°F), the readability 


~ of the recorder (£1°R)}, and the varintion in temperature 
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caused by inaccuracy in the placement of the thermocouples 


(+1°P, estimated). The overall uncertainty in the experi- 
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mental data was also calculated from the second power equa- 


? tion as 
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These uncertainty bounds are also shown in Figures 12 and 13. 
A procedure, similar to that used to find the uncertain- 
} 
| 
i 
! 
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WIRE * 
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ties of the analytical solution, was used to analyze the 
resulting uncertainty in the TRUMP numerical calculation, 


The results of these calculations are shown in Figures 14 
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and 15. fhe uncertainty in temperature varied with time 
. with a maximum variation of +2.95°F at the center of the 
rocket motor and a maximum variation of +1.95°F at the sur- 
? face of the motor. The time delay varied from +20 minutes 
at the center of the motor to +9 minutes at the surface of 
the motor. 
On the basis of the propagation of uncertainty analysis, 
it was determined that the solutions were most sensitive, in 
order of importance, to changes in the volumetric heat 


capacity, emissivity, and the conductivity. 


see oe 
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